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Abstract 

We study the variation of the broad-band spectral energy distribution (SED) of the BL Lac object MrkSOl as 
a function of source activity, from quiescent to flaring. Through -minimization we model eight simultaneous 
SED datasets with a one-zone Synchrotron-Self-Compton (SSC) model, and examine how model parameters 
vary with source activity. The emerging variability pattern of MrkSOl is complex, with the Compton com- 
ponent arising from 7-e scatterings that sometimes are (mostly) Thomson and sometimes (mostly) extreme 
Klein-Nishina. This can be seen from the variation of the Compton to synchrotron peak distance according to 
source state. The underlying electron spectra are faint/soft in quiescent states and bright/hard in flaring states. 
A comparison with Mrk421 suggests that the typical values of the SSC parameters are different in the two 
sources: however, in both jets the energy density is particle dominated in all states. 

Subject headings: BL Lacertae objects: general - BL Lacertae objects: individual (MrkSOl) - diffuse radiation 
- gamma rays: galaxies - 



1. INTRODUCTION 

The fueling of supermassive black holes, hosted in the cores 
of most galaxies, by infalling matter is thought to produce the 
spectacular activity observed in AGNs. In 10% of cases 
powerful coUimated jets shoot out in opposite directions at 
relativistic speeds. If a relativistic jet is viewed at a small 
angle to its axis, the observed emission is amplified by rel- 
ativistic beaming (Doppler boosting and aberration; see Urry 
& Padovani 1995). Sources, whose boosted jet emission dom- 
inates the observed emission (blazars0), are crucial to study- 
ing the physics of relativistic extragalactic jets. 

The jets' origin and nature are unclear They may be (ki- 
netic or electromagnetic) flows that dissipate some of their 
energy in moving regions associated with internal or external 
shocks. This complex picture is approximated, for the pur- 
pose of modelling the observed emission, with one (or more) 
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' Extreme blazars, whose thermal emission is intrinsically weak (i.e., no 
emission lines in their spectra), are called BLLac objects. 



relativistically-moving, magnetized, homogeneous plasma re- 
gion (blob), where a time-varying non-thermal population of 
electrons emit radiation (e.g., Maraschi et al. 1992). The lat- 
ter is synchrotron radiation and its comptonized (IC) coun- 
terpart, that peak at (respectively) IR/X-ray and GeV/TeV 
frequencies [in low/high-frequency -peaked BLLac sources 
(LBL/HBL)] and show correlated luminosity and spectral 
changes. In the case of HBLs one same population of elec- 
trons in the blob generates the low-energy (IR to X-ray) syn- 
chrotron electrons and Compton-(up)scatters them to high en- 
ergies [Synchrotron-Self-Compton (SSC) mechanism], with 
no external sources of soft photons. 

The recent availability of simultaneous broad-band SEDs 
(mostly for nearby BLLac objects) has enabled addressing 
the important issue of how the emission changes as a function 
of the source's global level of activity: i.e., given an emission 
model that fits the data, it is should be examined what model 
parameters are correlated with source activity. Another re- 
quirement concerns using a full-fledged -minimization pro- 
cedure instead of the "eyeball" fits still most commonly used 
in the literature, in order to obtain unbiased results. 

In a previous study of this type, Manzkuzhiyil et al. (201 1) 
examined nine SEDs of the HBL object Mrk421 (z=0.031). 
In this follow-up paper we analyze eight SEDs of the HBL 
source MrkSOl (^=0.034), applying the same one-zone SSC 
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emission model and -minimization procedure (Sect.2 and 
3.1, respectively) to the datasets described in Sect. 3. 2. The 
results are presented and discussed in Sect.4, and summarized 
in Sect.5. 

2. BLLACSSC EMISSION 

Following Mankuzhiyil et al. (2011), to describe the HBL 
broad-band emission we use a one-blob SSC model (Tavec- 
chio et al. 1998, later T98; Maraschi & Tavecchio 2003). This 
adequately describes broad-band SEDs of most HBLs (Tavec- 
chio et al. 2010) and, for a given source, both its ground and 
excited states (Tavecchio et al. 2001; Tagliaferri et al. 2008). 
Convincing evidence for the SSC model is the X-ray/VHE-7- 
ray variability correlation (e.g., Fossati et al. 2008): since in 
the SSC model the emission in the two bands is produced by 
the same electrons, a strict X-ray /7-ray correlation is expected 

0. 

In this work, for simplicity we use a one-zone SSC model, 
assuming that the entire SED is produced within a single ho- 
mogeneous region of the jet. Although generally adequate 
to reproduce HBL SEDs, one-zone models have troubles ex- 
plaining some specific features of blazar TeV emission. In 
particular, while very large Doppler factors are often required 
in a one-zone model, radio VLBI observations hardly detect 
superluminal motion at parsec scale (e.g.. Finer et al. 2010; 
Giroletti et al. 2006). This led to the proposal of the existence 
of a structured, inhomogeneous and decelerating emitting jet 
(Georganopoulos & Kazanas 2003; Ghisellini et al. 2005). 
Inhomogeneous (two-zone) models (e.g., Ghisellini & Tavec- 
chio 2008) have been also invoked to explain the ultra-rapid 
variability occasionally observed in TeV blazars (e.g., Aharo- 
nian et al. 2007; Albert et al. 2007). 

In the one-zone SSC model the emitting plasma is con- 
tained in a spherical blob of radius R in relativistic motion 
(described by a bulk Lorentz factor F) along the jet at an an- 
gle 6 w.rt. the line of sight to the observer, so that special rel- 
ativistic effects are cumulatively described by the relativistic 
Doppler factor, 6 = [F(l-/3cos6')]-i. The blob is filled with 
a homogeneous tangled magnetic field with intensity B and 
by a population of relativistic electrons of density rig, whose 
spectrum is described by a broken power-law (PL) function of 
the electron Lorentz factor 7, 



Neil) 



7min < 7 < 7br 
7br < 7 < 7ma3i 



(1) 



(Of course, rig = X^""'' ^e(7) d7 .) This purely phenomeno- 
logical choice is motivated by the observed shape of the 
humps in the SEDs (e.g., T98). 

In shaping the VHE 7-ray part of the spectrum it is impor- 
tant: (/) to use the full Klein-Nishina (K-N) cross section (see 
Maraschi & Tavecchio 2003) in the T98 model; and to 
correct > 50 GeV data for Extragalactic Background Light 
(EBL) absorption - as a function of photon energy and source 
distance (e.g., Mankuzhiyil et al. 2010): here we use the pop- 
ular Franceschini et al. (2008) EBL model. 

It is useful to sample the SED around both peaks, because 
the synchrotron emissivity contains degeneracy between B 
and Ke (i.e., oc KeB^'^+^^Z'^, with n = 2a + 1 where 
a n is the electron's spectral index and a is the measured syn- 

' The rarely occurring "orphan" TeV flares, that are not accompanied by 
variations in the X-ray band, may arise from small, low-B, high-density 
plasma blobs (Krawczynski et al. 2004). 



chrotron emission profile) that cannot be lifted without knowl- 
edge of the IC peak, which provides the desired second equa- 
tion (e.g., jic oc ii'g _B'"+^'/^). Therefore, in general only 
knowledge of observational quantities related to both SED 
humps enables reliable determination of all SSC parameters. 
In particular, the one-zone SSC model can only be fully con- 
strained by using simultaneous broad-band observations. Of 
the 9 free parameters of the SSC model, 6 specify the electron 
energy distribution (rie, 71, 7br, 72, ni, n2), and 3 describe 
the global properties of the emitting region (B, R, 5). Some 
observational quantities are linked to the SSC model parame- 
ters: the slopes, 0:1.2, of the synchrotron bump on either side 
of the peak are linked to ni_2 (whence the phenomenological 
choice of a broken-PL electron spectrum); the synchrotron 
and IC peak frequencies, 1^8,10, and luminosities, Ls.ic. 
linked with B, ric, 5, 7br; finally, the minimum variability 
timescale tvar provides an estimate of the source size through 

3. -MINIMIZATION 
3.1. The code 

We start by setting 71 = 1, as usually assumed in the 
literature. A search for the best values (and associated un- 
certainties) of the remaining 8 free parameters is performed 
by best-fitting our chosen SSC model to the SED data: the 
corresponding -minimization is based on the Levenberg- 
Marquardt method - an efficient standard for non-linear least- 
squares minimization that dynamically interpolates between 
two different minimization approaches, the inverse Hessian 
method and the steepest descent method. 

A description of how the code works is given in Manku- 
zhiyil et al. (2011) and Ansoldi (2012). Let's stress two 
important points in our implementation. ( /) The T98 model 
only returns the numerical sample of the SSC SED. On the 
other hand, at each step of the loop the calculation of re- 
quires evaluating the SED at all the observed frequencies. If 
the model function is known analytically, these evaluations 
are a straightforward algebraic process. In our case, instead, 
we know the model function only through a numerical sam- 
ple coming from the implementation of the T98 model, and 
in general these sampled points do not include the observed 
points. The latter, however, typically fall between two sam- 
pled points, which - given the dense sampling of the model 
SED - allows us to use interpolation to approximate the SED 
value at the observed point with negligible uncertainties. 
The Levenberg-Marquardt method requires the calculation of 
the partial derivatives of with respect to the 8 floating SSC 
parameters. Contrary to the case when knowledge of the an- 
alytical form of the model function enables all derivatives to 
be obtained in closed form, in our case the derivatives are nu- 
merically obtained by evaluating the increment of with re- 
spect to a sufficiently small, dynamically adjusted increment 
of each parameter As this procedure could have led to some 
computational inefficiency due to the recurrent need to eval- 
uate the SED from the model at many finely-spaced points in 
parameter space (a CPU-time consuming operation), we set 
up a specific algorithm to minimize the number of calls to 
T98 across different iterations. 

Because of the subtleties in non-linear fits we optimized the 
procedure to identify the parameters uncertainties for reliabil- 
ity and statistical significance. The standard approach (uncer- 
tainties as the square-root of the covariance matrix diagonal 
elements) depends on a reUable quadratic approximation to 
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Figure 1. Eight SEDs of Mrk501. See Table 1 for details. 



the -surface around the minimum: otherwise uncertainties 
can be underestimated/overestimated. To avoid this we di- 
rectly solve for them using (we call the SSC parameters pi 
here) 



X^(pi"""\ . . . . . . ,P8""'0 =Xmin + "7 i= 1,2,...,8, 

where a = 1, 2.71, 6.61 reahze a 68%, 90% and 99% confi- 
dence level, respectively^ and x^in — X^(Pi™'"\ • • ■ jPs™™^) 
is the chi-square at the optimal estimation for the free pa- 
rameters. Starting from the minimum point, the above 
equations are solved twice for each parameter pi, keep- 
ing the other seven fixed to their optimal values. We 
first add to the floating parameter positive increments. 
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Pl ' and p\ represent what is called a bracketing of the 
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solution pl ' at which x \V\ i - ■ ■ iVX t ■ ■ iPs ) = 
Xmin + 2.71. Given the bracketing, we numerically solve for 
p[°'''^ using Ridders' method (Press 1992) and define the 
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and the left uncertainty is then al ' ~ Pl ' — Pl 
The final results for the parameters are then estimated as 

fmin)+<T'"^''" 

Pi (left) the 90% confidence level. 

' Our choice is the standard 90% confidence level. 

* Ridders' method (Ridders 1979) only requires evaluation of the func- 
tion which is searched for a zero, not of the function deiivatives: it is thus 
well suited to our case, in which calculation of (and thus model function) 
derivatives is computationally demanding. The order of Ridders' method is 
\/2 (Press 1992), so it is superlinear and competitive with more refined ap- 
proaches. 



3.2. The data sets 

Searching the literature for simultaneous broad-band SED 
data sets of MrkSOl returned the following: 

• State 1 data were collected by Suzaku (March 23-25, 2009), 
by MAGIC and VERITAS (March 23 and March 24 2009, re- 
spectively), and by Ferm//LAT (all throughout the campaign) 
- see Acciari et al. (2010). 

• State 2 data were taken from the most intensive multi- 
wavelength campaign ever conducted on Mrk501 (Abdo et 
al. 2011). Swift, RXTE, Femii/LKT, MAGIC and VERITAS 
data are used in this work. The X-ray data were obtained dur- 
ing Marl5-Aug 1, 2009. Fermi observed the source in the 
same time window, but data taken in May were omitted in 
this work due to a variable high-energy flux. MAGIC observa- 
tions were carried out between Apr27-May 9, 2009: however, 
due to bad weather and a shutdown for hardware upgrade, the 
observation days were significantly reduced from the orig- 
inal schedule. VERITAS observations took place between 
Marl7-Jun22, 2009: but data from May 1-4 were dropped 
due to high variability. 

• State 3 data come from a multi-frequency campaign in July 
2006 (Anderhub et al. 2009). KVA (optical) and MAGIC 
(VHE7-ray) observations were carried out on July 18, 19, 
and 20 - plus a few tens of minutes in August and Septem- 
ber 2006. Suzaku observations were carried out on July 18 
and 19, 2006. During the whole campaign, no significant flux 
variability was observed in any of the three frequency ranges. 

• State 4 and State 5 data come from a campaign carried out 
by RXTE and HEGRA in June 1998: State 4 uses observa- 
tions from June 15-26, 1998, whereas State 5 uses data taken 
during a higher-activity state on June 27-28, 1998 (Sambruna 
et al. 2000). 

• States 6, 7, 8 are based on a multi-frequency campaign, in- 
volving CAT and BeppoSax, from Mar-Oct 1997. The X-ray 
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data of states 6 and 7 were obtained on April 7 and 11, re- 
spectively. The VHE 7-ray fluxes used in State 6 and 7 rep- 
resent average emissions of the low and intermediate states 
observed during the campaign (flux < 1.2 x 10^^° cm~^ s~^ 
and 1.2 X 10"™ cm^^ g-i < flux < 5.3 x 10"™ cm^^ s-\ 
respectively). State 8 comprises X-ray and VHE 7-ray data 
obtained during the giant flare of April 16, 1997 (Djannati- 
Ataietal. 1999). 

4. RESULTS 

Best-fit SSC models of simultaneous blazar SEDs are cru- 
cial to measuring the SSC parameters that describe the emit- 
ting region. The emerging best model is obtained by thor- 
oughly searching for the absolute minimum. Our proce- 
dure ensures there is no obvious bias affecting the resulting 
best-fit SSC model. 

The SED data and their best-fit SSC fit parameters are re- 
ported in Table 1 . The corresponding fits are plotted in Fig. |2] 
The 90% error bars reported in Table 1 result from a detailed 
study of the profiles around the minima corresponding to 
each parameter (see Sect. 3. 2), and in general they turn out 
asymmetric. Within the available statistics (8 data points), 
source activity (measured as the total luminosity of the best- 
fit SSC model0 seems to be most correlated with the electron 
spectrum normalization, (see Fig. O, and with the high- 
energy power-slope, n2' both quantities vary with luminosity, 
in the sense that an increase of luminosity means an increase 
of the electron density and its high-energy fraction. 

In more detail, from Table 1 - within the low statistics in- 
volved - we can see that: 

(i) in quiescence (states 1, 2, 3) apparently uncorrected 
variations of the SSC parameters produce SEDs of constant 
(isotropic) luminosity; 

(ii) in the sequence of June 1998 (states 4, 5) a '-~^3-fold in- 
crease of luminosity is accompanied by a 3-fold increase of 
Ke and by a substantial electron spectral hardening - as well 
as by a 2-fold decrease of 6; 

(Hi) in the Giant Flare of 1997 (states 6, 7, 8) the 4-fold in- 
crease of luminosity is accompanied by increments of and 
7br by factors of, respectively, 5 and 2, and by a strong hard- 
ening of the electron spectrum - as well as by a decrease 
of S that eventually amounts to a ~2-fold decrease in rela- 
tivistic boosting. Along this sequence the emission radius R 
increases, and the magnetic field B decreases, such that the 
quantity BR^ - i.e., the magnetic flux across the boundary of 
the emitting region - appears to be approximately conserved; 

(iv) although in most emission states considered in this paper 
n2 — ni ^ 0, in states 8 and 5, i.e. in the flaring culminations 
of the April 1997 and June 1998 periods of activity, n2 has 
decreased such that 712 — ni ^ 0.5: i.e., the electron spec- 
trum, generally assumed to be double-PL, tends to look like a 
single-PL spectrum during flares. 

' The (isotropic) bolometric luminosities used in tliese plots have been 
obtained directly from the model SED. After determining the numerical ap- 
proximation to the SED log[uF{u)], with the parameters being fixed at their 
best values obtained with the previously described minimization procedure, 
we have peiformed L = 47rD^ J^'"'"' F{v)du, with the luminosity 
distance, and !^niin. '^max set at 2.5 and 0.75 decades, respectively, below the 
synchrotron peak and above the Compton peak. In this way we make sure 
to perform the integral over all the relevant frequencies in a way that is inde- 
pendent from the location of these peaks. (This procedure was also used in 
Mankuzhiyil et al. 20 11 to calculate the luminosity of Mrk421, in spite of a 
typo that makes the description inaccurate.) 



In general, we find that the kinetic energy density of rel- 
ativistic electrons, Ue = ni^c? J^'"'"' j Ne{j) dj which is a 
combination of K^, fii, n2, 7min, Jhr, and 7niax, is strongly 
correlated with the blob rest-frame luminosity (Fig. |3)- 

4.1. Electron injection spectrum 

The large variation in luminosity and SED shape displayed 
by Mrk501 between flaring and non-flaring states is matched 
by a concomitant change in the electron spectrum, N^i-f) (see 
Fig.©: 

(i) In the three quiescent states (Fig. [3] left) the spectra are 
clearly double-PL (An 2± 1.6) with similar 7br ~ 6 x 10^ 
and normalization ranging from ~70cm"^ to 230 cm~'^ 
to 255cm"3. 

Along the sequence of the two states of 1998 June 15-26 
and 27-28, the spectrum hardens (An from 1.32 to 0.61) and 
brightens (Ke increases by a factor of 3) (Fig.|3] center). 
(Hi) Along the sequence of the three states (6,7,8) from 1997 
April 7, 11, 16 the spectra become harder in shape and higher 
in normalization as the luminosity rises (Fig. m right). On 
April 7 the normalization is lowest (Ke — 90 cm"'^ and the 
spectral break is highest (An = n2 — ni — 1.35). By April 11 
the normalization is intermediate (Ke — 165 cm~^) and the 
spectral break has decreased (An = 1.10). On April 16, at 
the peak of the flare, the normalization is highest (Ke — 465 
cm"^) and the spectral break has almost disappeared (An = 
0.38). 

The implied temporal pattern is quite clear At the onset 
of the flare, the electron spectrum has a distinctly double-PL 
slope and relatively low normalization. As the flare goes on, 
more and more newly accelerated electrons are injected, and 
the resulting spectrum is a combination of the "old", double- 
PL spectrum and the "new", single-PL injection spectrum. 
At the peak of the flare the injection spectrum dominates the 
emission. Correspondingly, the electron energy density (Ue) 
shows an increasing trend with luminosity and, within a cor- 
related sequence of time-ordered states, with the development 
of the flare (see Fig. |3]l. 

4.2. Thomson versus extreme K-N scattering regimes 

We can check whether the elementary scattering processes 
that give rise to the IC component of the SED occur mainly 
in the Thomson or extreme K-N regime by studying the log- 
arithmic distance between the synchrotron and IC peaks. At 
its simplest, the argument goes as follows. In the electron's 
rest frame (primed quantities), by scattering off the electron a 
photon with initial energy e' goes off at a scattering angle 9[ 
and an energy e[ given by 



^ l + {e'/meC^){l-cose[) ' ' 

(e.g., Blumenthal & Gould 1970). In the lab system the en- 
ergy is 

ei = 7e'i[l + /3cos(7r-6i'i)] = -fe[ (1 - (3cos9[) , (^) 

so ei ^ 7ep Therefore, from Eq.(l2} we see the two limiting 
cases: 

(/) in the low-frequency limit e' ^ TOgC^ (Thomson regime), 
it is e[ e' . The energy transformation from the lab frame to 
the electron frame, e' = 76 (1 — /3cos0), implies e' ~ je, so 
d ~ 7^e. In the SSC framework ei = hi'ic and e = hi^s, and 
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Figure 2. Variations of the plasma blob's radius (R), magnetic field (B), and Doppler factor (5), and of the electrons' spectral indices (ni, n2) and spectral break 
(7br)- Also shown are the variations of the magnetic and electron energy density (Ub and Ue, respectively). Circles and triangles denote Mrk501 and Mrk421, 
respectively. Uncertainties for Mrk421 (Mankuzhiyil et al. 2011) have been estimated again with the method described in this paper to make them consistent 
with those of Mrk501. 



the low-frequency condition writes 7 ^ rrieC^ / [hvs); since 
z/s = 5 X 10^725 Hz (in the blob's frame)[3and B ~ 0.02 G 
in all the eight Mrk501 states under consideration, the latter 
condition translates into log 7 <C 5. We then have hence 



log(;^ic/i^s) 
log7 



log7 < 5.0 ; 



(4) 



(ii) in the high-frequency limit e' ^ mgC^ (extreme K- 
N regime) it is e'^ = rUeC^, so ci ~ jnieC^- In the 
SSC framework we have, for one electron, ei = /ific, so 
Vie - iruec^h-^, and = 2.8 x 10^7^5 02. For the 
case, of relevance here, of a broken-PL electron spectrum, it 
is Vie — l'meC^h^^g{ai-,oi2) [with f;(ai, 02) — exp[(ai — 
l)-i+0.5(a2-ai)"^], being ai^2 = (ni,2-l)/2] (seeT98) 
and Vs = 'ox 10^725 Hz. For values B = 0.02 G, ni = 1.7, 
71.2 = 3.0, appropriate for MrkSOl (see Table 1), we get 



\og{viclvs) 
log7 



14.75 

log7 



1 , log7 » 5.0 ; (5) 



The predictions in Eqs.(4) and (5) can be graphically tested 
by plotting (in a log-log plane) the IC-to-synchrotron peak 

For a single electron it is Us = 2.8 X (e.g., Tucker 

1975). For the electron distribution in Eq.llT) the average value of 7 is 
7 = f{ni, n2) 7br. Typical values of ni, n2 imply 7 ~ 1.37br- 



distance vs. 7 (see Fig. I4.2in states 1-5 (quiescent to mod- 
erately active) are located on, or close to, the Thomson line 
whereas states 6-8 (active) are located close to the extreme 
K-N curve. 

4.3. Parameter and source evolution 

Examining Table 1 , we may tentatively point out some 
trends shown by the SSC parameters during the evolution of 
an active event. States 6, 7, and 8, snapshots of the active 
event of April 1997 that culminated in the Giant Flare (Pian 
et al. 1998), taken At days apart from one another, suggest 
that the quantity B is conserved (as expected in a fully 
ionized plasma), and - assuming the rise to the Giant Flare 
to have started At days before state 6 - that the blob radius 
grows with time as i? cx t^/^ . 

If real, the latter suggested regularity corresponds (Wand & 
Kusunose 2002) to the self-similar solution (Sedov 1959) of 
the subrelativistic expansion of a plasma blob. Pursuing this 
analogy further, this Sedov expansion leads to the formation 
of a strong shock which then sweeps some electrons from the 
medium and accelerates them to relativistic energies: so the 
expansion acts as the injection of relativistic electrons (Wand 

' ' Since in our assumed scheme the dominant synchrotron power is emit- 
ted by electrons at the break (i.e., with 7br): and the dominant IC power (in 
the Thomson regime) is accounted for by photons with frequency ~ I'iu-^s, 
the relevant 7 in this argument is 71-,^ ■ 
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Figures. Electron spectra (top) and energy densities (bottom) in Mrk501's jet during: (left) quiescent times (uncorrelated states 1, 2, 3); (middle) June 
1998 event (time-ordered states 4, 5); (right) April 1997 event (time-ordered states 6, 7, 8). The electron spectra, N(j), are plotted here multiplied by 7^ in 
order to better illustrate elecrons of which energy carry the bulk of the kinetic energy. In the assumed one-blob SSC model, the quantity L<5~* represents the 
bolometric SSC model luminosity corrected for the effects of relativistic beaming (see Urry & Padovani 1995): i.e., L<5~* represents the blob-frame isotropic 
SSC luminosity. 



& Kusunose 2002). 

These elements may suggest the following picture for the 
early development of a flare: the blob expands while traveling 
along the jet, and fills it because its internal radiation pressure 
swells it until magnetic confinement forces prevail and push 
the blob further out where it can keep expanding. As long 
as new particles are accelerated so that the internal pressure 
builds up, one can conjecture a positive feedback among ex- 
pansion, motion along the jet, internal shock formation, and 
particle acceleration until maximum emission is reached. 

4.4. Comparisonwith Mrk421 

A cursory comparison with Mrk421 (Mankuzhiyil et al. 
2011) suggests that MrkSOl (see Table 1) has a larger (by a 
factor ^12) emission size R; a lower but more stable Doppler 
factor ((5 - 20 ± 5) [whereas Mrk421 has <5 - (25 ± 5) - 
(100 ± 20)]; and a lower B. So the two jets appear to be 
substantially different from each other. 

A deeper insight into similarities and differences in the 
emission physics of the two jets can be reached by examin- 
ing congelations between SSC parameters (e.g., Tavecchio et 
al. 1998; see Fig. |4]l. A i3-7br anticorrelation is predicted 
if the synchrotron peak, oc B^^, remains (roughly) con- 
stant from state to state. For constant and v\c, a B-6 re- 
lation is predicted to be inverse in the Thomson regime and 
direct in the extreme K-N Hmit (e.g., Tavecchio et al. 1998). 



Both B-^j\,i and B-5 inverse correlations do hold for the 9 
SEDs of Mrk421 examined by Mankuzhiyil et al. (201 1). For 
MrkSOl the frequencies and v\c are not constant, hence the 
B-'fhi (anti)correlation does not hold (see Fig. |4]i. Similarly, 
in the B-S plot there seems to be no correlation for MrkSOl, 
possibly suggesting a mix of (predominantly) Thomson and 
(predominantly) extr eme K-N states for this source. More di- 
rectly, from Fig. l4.2l we see that all Mrk421 points are on, or 
close to, the Thomson line, where the synchrotron and Comp- 
ton peak frequencies are related by i^c l^s (x 7^, whereas the 
MrkSOl points are found close to both the Thomson and the 
extreme K-N curves. 

Using the best-fit SSC parameters, for each SED of both 
galaxies we compute the electron energy density C/e and the 
magnetic energy density Ub- The results (see Fig. |4]i show 
that the the situation is far from particles-field equilibrium, 
with the electrons dominating over the field by orders of mag- 
nitude. □ (If we would also consider a contribution from 
accelerated protons, which appreciably contribute not to the 
radiative yield but to the mechanical power of the jet, our con- 
clusion is even stronger.) These results suggest that in all the 

In both sources C/e is computed assuming 7niin = 1- The impact of 
7jnin on Ve is not dramatic as long as n\ < 2, which is the case here. So our 
conclusion on particle dominance, although based on one particular value of 
7niin: retains general validity. 
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Figure 4. Logarithmic distance between the IC and synchrotron peaks as a 
function of the electrons' spectral break energy for the 8 SEDs of MrkSOl 
(filled circles) considered in this paper and the 9 SEDs of Mrk421 stud- 
ied in Mankuzhiyil et al. (2011; empty squares). The thick solid and 
dashed lines correspond to the Thomson (hu/rrieC^ — 0) and extreme K- 
N (hu/meC^ — oo) limits, respectively. In calculating the extreme K-N 
curves magnetic field values of B = 0.02 G (appropriate for MrkSOl, see 
Table 1), and B = 0.4 G and 1 G (appropiiate for Mrk421, see Mankuzhiyil 
et al. 2011) were assumed. (In all cases, when not indicated otherwise, we 
used «! = 0.5, 02 = 1.5.) Thin lines represent the extrapolation of the two 
limiting regimes into the hu ^ rrte region. In the extreme K-N limit, the 
two peak frequencies are closer than in the Thomson limit. 



States considered in this paper, quiescent and active, the jets 
of both Markarians seems to be primarily kinematic - as op- 
posed to electromagnetic. This result was already pointed out, 
for quiescent states alone, by e.g. Kino et al. (2002). Both 
the magnetic and the electron energy densities are distinctly 
higher in Mrk421 than in MrkSOl; however Mrk421 sug- 
gests a Ub oc C/g behavior, whereas no correlation appears 
for MrkSOl - especially without the two lowest- [/ b, highest- 
Ue points that correspond to the flaring states S and 8 (see 
Fig.H. 

5. SUMMARY 

Through -minimization we have modeled eight simulta- 
neous SED datasets of the HBL source MrkSOl with a one- 
zone SSC model. Our results suggest that MrkSOl's vari- 
ability is complex, with the IC component being produced 
by photon-electron scatterings that at times are in preva- 
lent Thomson regime and at times in prevalent extreme K-N 
regime. The variations of the SED have enabled us to track 



the corresponding variations of the underlying emitting elec- 
tron spectrum: from distinctly broken-PL and less intense and 
softer in the non-flaring states, to virtually single-PL and more 
intense and harder at the peak of flares. The former are inter- 
preted as the steady-state spectra of aging electrons, the latter 
as the fresh electron injection spectrum. Comparing data for 
MrkSOl and Mrk421, the two jets mainly differ by the phys- 
ical characteristics of their emission regions, but their energy 
budgets appear to be both particle dominated. 
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Table 1 

Datasets and best-fit single-zone SSC model parameters for the eight datasets of MrkSOl. States are named as in Fig.l. References are: [1] Acciari et al. (2010), 

[2] Abdoetal. (2011), [3] Anderhub et al. (2009), [4] Sambruna et al. (2000), [5] Djannati-Atai' et al. (1999). B is given in units of 10" ^G, ii in units of 
IQi^cm, Kin units of cm~3^ ^y^^ in units of 10*, 7max in units of 10^. The isotropic luminosities (in erg s~i, given in log), and the synchrotron and IC peak 

frequencies (in Hz, given in log), are derived from the SEDs' best-fit SSC models. 



